Helical Magnetic Fields from Sphaleron Decay and Baryogenesis 
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Many models of baryogenesis rely on anomalous particle physics processes to give baryon number 
violation. By numerically evolving the electroweak equations on a lattice, we show that baryogen- 
esis in these models creates helical cosmic magnetic fields. After a transitory period, electroweak 
dynamics is found to conserve the Chern-Simons number and the total electromagnetic helicity. We 
argue that baryogenesis could lead to magnetic fields of nano-Gauss strength today on astrophysical 
length scales. In addition to being astrophysically relevant, such helical magnetic fields can provide 
an independent probe of baryogenesis and CP violation in particle physics. 



Cosmic magnetic fields can arise from a number of 
early universe processes {e.g. P, 0, |3, 0, H, The 
detection of such magnetic fields would be an important 
step in our understanding of structure evolution, and the 
magnetic field distribution could be used as a probe of 
early universe cosmology. Most magnetic field generation 
mechanisms discussed so far, produce non-helical fields, 
and thus it is of particular interest that processes like 
electroweak baryogenesis imply magnetic fields that are 
helical @|. If helical primordial magnetic fields are 
observed, they would offer confirmation of cosmological 
baryogenesis, and the magnetic field properties could be 
turned into a detailed probe of particle physics and cos- 
mology at the epoch of baryogenesis. In particular, the 
CP violation that leads to a universe filled with matter 
and no antimatter would be probed by cosmological ob- 
servations. 

These considerations have prompted us to study 
the production of magnetic fields during the decay of 
sphalerons 0, a process that is key to electroweak 
baryogenesis. While electroweak baryogenesis is not 
successful in the minimal standard model, our study 
also applies to any extension of the standard model in 
which baryon number violation is mediated by sphaleron- 
like processes involving changes in the winding (Chern- 
Simons number) of vacuum non-Abelian gauge field con- 
figurations. 

To understand the connection between sphalerons and 
helical magnetic fields, it is simpler to think of "deformed 
sphalerons" where the gauge field configuration resembles 
that of electroweak strings The sphaleron can then 
be interpreted as linked loops of electroweak Z-string 
or a confined electroweak monopole-antimonopole pair 
[13, E Q El- The linked loops of Z magnetic flux can 
then decay into linked electromagnetic flux as described 
in and thus the resulting electromagnetic field carries 
magnetic helicity. If we think of the sphaleron in terms 
of the confined magnetic monopole pair, the string that 
confines them is twisted and this also leads to magnetic 
helicity. In [^, 0| these considerations indicated a remark- 



ably simple relationship between the cosmic magnetic he- 
licity density and the cosmic baryon number density 
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where we consider some large spatial volume V, A is the 
electromagnetic vector potential, is the average baryon 
number density and a — 1/137 is the fine structure con- 
stant. 

Our goal is to examine the heuristic arguments in 
[E Hi by explicitly studying the decay of an electroweak 
sphaleron. (Recently, along similar lines, magnetic fields 
produced during preheating at the electroweak scale were 
studied in Ref. [7[.) We will indeed find that helical 
magnetic fields are generated during sphaleron decay and 
the relation in Eq. ([T]) holds at the order of magnitude 
level. Our numerical results also show, somewhat unex- 
pectedly, that the Chern-Simons number and the electro- 
magnetic helicity are conserved after a transitory initial 
period (also see [1]). We also reconsider the net mag- 
netic field generated during baryogenesis and find that 
the field strength is likely to be much larger than has 
been previously estimated. 

We work with the bosonic sector of the electroweak 
Lagrangian 



where 
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D^^ = [d^ - I'-B^ - V^T ■ W^j $, (3) 

The SU{2) generators, r", are normalized by Tr(r°r^) = 
2(5"''. The U{1) (hypercharge) field strength is 

B^, = d^B, - d,B^, (4) 

and the SU{2) field strength is defined as 

W^;. = d,W^ - a.iy; -f ge'^'^W'^w: (5) 
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To study the magnetic field produced by a decaying 
sphaleron, we numerically set up a configuration like the 
electroweak sphaleron solution [l^, a per- 

turbed sphaleron, since the exact solution, even though 
it is unstable, will take a long time to decay numerically. 
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where 



Ge {x) = exp [ie(r)T • x/2] 



(6) 
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is an SU{2) gauge transformation, with 8(r) = 
TT (1 — exp(— r)), ensuring that the gauge fields fall off 
fast enough away from the origin. The functions f{£,), 
/i(^) and p(^) are profile functions and are taken to de- 
pend on the radial coordinate ^ — gvr j \f2 alone. 

In the true solution, the profile functions also depend 
on the angular coordinates and have to be determined 
numerically [17| . Since we want to start with a perturbed 
sphaleron we do not need the precise forms of / and 
/i, but we do need to meet the asymptotic properties 
satisfied by / and h: /, ft, ^ as r ^ 0, and f,h^ 
1 as r — > CXI. Following Klinkhamer and Manton [lo| . 
we choose for the profile functions their Ansatz, labeled 
a, with the length scales given by S = 3.79 and fl = 
1.90. We take for the initial U{1) gauge field, Eq. ([8]), 
the small g' approximation [lol |. The profile function p 
can be found, in terms of / and h, by solving the U{1) 
equation of motion in the SU (2) background. 

The decay is studied by evolving the SU{2) x [/(I) 
electroweak field equations using the standard Wilsonian 
approach for lattice gauge fields 18, 19, 20, 21, 22]. The 
temporal gauge, = Bo = 0, allows a simple identifica- 
tion of the canonical momentum, and we use it through- 
out in this work. We adopt Graham's implementation of 
the lattice equations [23 | 



(l){t+) = 2(l>{t) - (f>{t~) + (Atf 
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where the gauge fields are on the links and are defined in 
terms of the matrices by 



jjP ^ ^ig'BJAx/2^igWf-rAx/2 



(10) 



Here p labels the point in the lattice and j the link em- 
anating from point p. The U matrix along a timelike 



link, j = t, is equal to the identity matrix in the tem- 
poral gauge. For the link from site p in the negative j 
direction, we take: 

The evolution of the gauge fields is given by 



= exp 
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logC/g 



iAt^Ax {g'jP + gjP 



(12) 



Here, we define U"^--, 



UPUP+ojjP_+3+fjjP_+f fo^ the 



plaquette and JJ = -g'luv{(j>P'^U^(t>P+^) /Ax and JJ = 
— g'lm (</)^''l'TC/j'0P+-') I Ax for the gauge currents. For 
further details, the reader is referred to [13] ■ In addi- 
tion to these evolution equations we have implemented 
absorbing boundary conditions as described in [2^ but 
extended to non-Abelian fields. 

The evolution of Chern-Simons number is of particular 
interest to us since it is correlated with changes in the 
baryon number. It is given by 
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where Np = 3 is the number of families and there is no 
surface term because our fields vanish at infinity. 

In the electroweak model, the electromagnetic gauge 
field is defined in terms of the electroweak W^^ and 
B'^ gauge fields 



A^ = svae^n^Wj^ + cos 
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(14) 



with n° = — ^V^^/^ttj). However, there is a choice 
of definitions for the electromagnetic field strength. For 
example. 



A^^ = sin^ 



COSf 



B„ 



(15) 



is the natural definition to calculate the energy density in 
the massless electromagnetic field. Yet this definition is 
not simply the curl of the gauge field in Eq. because 
derivatives of also involve derivatives of n". Even at 
late times, when becomes small, such derivatives 

will in general contribute. The only clean resolution of 
this issue that we have found is to define the electromag- 
netic gauge field in unitary gauge and then define the 
electromagnetic field strength as the curl of the gauge 
potential as in the usual Maxwell theory. At early times. 
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when Df^^ is significant, V x A does not coincide with 
Aij in Eq. ITSl) . but they do coincide at late times. Then 
the hehcity tells us something about the topology of the 
very same magnetic field lines that carry energy density. 

To go into unitary gauge, we calculate Af^ at every time 
step only after applying the SU(2) gauge rotation 
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where $i and $2 are the upper and lower components of 
$. It is easily verified that 52$ = V$t$(0, l)"^, det52 = 1 
and gl = 5^^. After the rotation, n° = (0,0,1), and 
when the rotation is applied to Wj^ the inhomogeneous 
part of the gauge transformation contributes in a non- 
trivial way to Af^. At late times, we then find that 
e^^'^Ajk/2 coincides with (V x A)*, and we can unam- 
biguously keep track of the total electromagnetic helicity 
defined as Hit) = Vh(t) (see Eq. ([J)). 

In Fig. [T] we show the time evolution of the Chern- 
Simons number and the total electromagnetic helicity. 
We see changes in the Chern-Simons number as well as 
growth in the magnitude of the electromagnetic helicity, 
and then a period of conservation (time step 50 to 250). 
While electromagnetic helicity is known to be conserved 
in a conducting plasma, its conservation in the present 
situation is novel because we are solving vacuum equa- 
tions and the only charges in the system are due to the 
fields. 

The value of electromagnetic helicity (w 2.5) is much 
less than the value estimated (w 200) in Refs. [1,0]. The 
estimate in [6|] assumed a certain decay channel for linked 
loops of Z-string and it is likely that the electromagnetic 
helicity is a function of the precise instability through 
which the sphaleron decays. A more rigorous estimate 
of the final helicity needs to be done statistically, taking 
into account the conditions at the epoch of baryogenesis. 

Once electromagnetic field helicity is produced, it will 
evolve, and eventually get frozen-in in the highly con- 
ducting ambient plasma. This is the scenario envisaged 
in Ref. Q and, under certain assumptions about the in- 
verse cascade, leads to the estimate that the cosmic mag- 
netic field is ~ 10~^^ G at recombination and coherent 
on a comoving scale of ~ 0.1 pc. We now argue that this 
is really an underestimate of the magnetic field strength 
and typically we can expect a much higher field strength. 

The essential point is that every baryon number vio- 
lating reaction goes via the sphaleron and produces mag- 
netic fields, whereas the estimate in only accounts for 
the magnetic field produced due to the excess of baryons 
over antibaryons that we see today. To make this clearer, 
suppose that sphalerons decay in some volume to produce 
Nt baryons while others decay to produce Nb — — e 
antibaryons, where e is entirely due to fundamental CP 
violation. Magnetic fields will be produced in each one of 
these Nb + Nb = 2Nb — e sphaleron decays. The baryon 
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FIG. 1: Time evolution of Chern-Simons number (dashed 
curve) and electromagnetic helicity (solid curve). The sign of 
the Chern-Simons number has been flipped for greater clarity. 



excess, however, is due to CP violation and is given by 
Nb — N\j = e sphaleron decays. The magnetic field pro- 
duced by just these excess number of reactions is much 
smaller. However, just as the baryons and antibaryons 
can annihilate, it is likely that some of the magnetic fields 
produced due to baryon and antibaryon production will 
also annihilate. It is hard to estimate exactly how much 
magnetic field survives, but the estimate in [6| is the min- 
imum value of the magnetic field. This value is protected 
by helicity conservation. 

We will now obtain another estimate with less restric- 
tive, and more realistic, assumptions about the evolution 
of the magnetic fields. Let us suppose that the mag- 
netic fields due to baryon and antibaryon production do 
not completely annihilate, leading to a magnetic field en- 
hancement by a factor r. Then 
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Successful big bang nucleosynthesis (BBN) constrains the 
energy density in the magnetic field so that B{t^cc) < 1 
G. We choose the BBN bound as opposed to other bounds 
since BBN places strong constraints on magnetic fields 
with small coherence scale e. 



see 
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(Other con- 
straints due to gravitational wave production may apply 
to the non- helical component of the magnetic field [25|.) 
This implies that 



r < 10 
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(18) 



Now to connect to particle physics, we realize that r de- 
pends on the total number of sphaleron processes, which 
are Nb + Nb in number. The estimate in of lO^^'^ 
G at recombination assumes that the magnetic field is 
proportional only to the net baryon number Nb — Nb- 
Assuming that the magnetic field strength gets reduced 
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due to annihilations in some stochastic way, we expect 



Nb-Nb 



(19) 



where 7 is an undetermined exponent which is 1 for no 
annihilation and 1/2 for Brownian evolution of the mag- 
netic field strength. 

The baryon excess is purely due to CP violation in the 
particle physics responsible for baryogenesis. For exam- 
ple, in the case of the electroweak model [26.] . 



Nb + Nb 
Nb~Nb 
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leading to 10^^ G at recombination for the stochastic 
case (7 — 1/2) and 10^ G for the no annihilation case 
(7 = 1). The latter estimate clearly violates the BBN 
bound in Eq. (fT8|) . while the former leads to 10~^ G 
today though on short coherence scales (0.1 pc). 

The argument above is intended to show that it may 
be possible to derive important model-independent con- 
straints on particle physics from limits on cosmic mag- 
netic fields. Further, if cosmic magnetic fields are ob- 
served, they can be used to derive detailed information 
about processes at baryogenesis and hence about high 
energy particle physics and CP violation. 
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